Melting of a 2D Quantum Electron Solid in High Magnetic Field 
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The melting temperature (T m ) of a solid is gen- 
erally determined by the pressure applied to it, 
or indirectly by its density (n) through the equa- 
tion of state. This remains true even for helium 
solids 1 , where quantum effects often lead to un- 
usual properties 2 . In this letter we present exper- 
imental evidence to show that for a two dimen- 
sional (2D) solid formed by electrons in a semi- 
conductor sample under a strong perpendicular 
magnetic field 3 (B), the T m is not controlled by 
n, but effectively by the quantum correlation be- 
tween the electrons through the Landau level fill- 
ing factor v=nh J eB . Such melting behavior, dif- 
ferent from that of all other known solids (includ- 
ing a classical 2D electron solid at zero magnetic 
field 4 ), attests to the quantum nature of the mag- 
netic field induced electron solid. Moreover, we 
found the T m to increase with the strength of the 
sample-dependent disorder that pins the electron 
solid. 

Electrons are expected to crystalize into a solid (so 
called "Wigner crystal"— ) when the (Coulomb) interac- 
tion energy between the electrons sufficiently dominates 
over the kinetic energy. One example of such an electron 
solid was found in a very low density (n) two dimensional 
electron system (2DES) realized on helium surfaces^ (at 
zero magnetic field). Because of the low n, the zero- 
point motion (given by the Fermi energy Ef = nh 2 /2irm 
where m is the electron mass) is negligibly small and 
at finite temperatures (T), as in the experiment^, the 
kinetic energy originates mainly from the classical ther- 
mal motion (fe#T). The melting of such a "classical" 
2D electron solid is determined only by the competition 
between the thermal kinetic energy and Coulomb inter- 
action (e 2 y / n/47re) [MKS unit is used exclusively in this 
paper] and is thought to be describable by the Kosterlitz- 
Thouless theory of 2D melting 6 . Experimentally, the 
melting was found to occur— at T cm = e 2 y/n/(AnekBT) 
(where e is the dielectric constant and fee the Boltzman 
constant) with r^l30, in excellent agreement with the- 
oretical calculations^. 

The 2DES as realized in high-quality GaAs/AlGaAs 
structures in our experiment has relatively high n, thus 
(in the absence of magnetic fields) the zero-point mo- 
tion (Ef) is significant and the 2DES does not solidify 
even at T—0. However, it is well known that a suf- 
ficiently strong perpendicular magnetic field (B) tends 
to suppress the kinetic energy of 2D electrons and in- 



duce the solidificatior&i .. On the other hand, at finite 
B the motion of electrons is quantized into Landau levels 
(LL) and delicate many-body quantum correlations^ 1 * 1 ^ 
among electrons can cause the 2DES to condense into 
fractional quantum Hall (FQH) 13 liquid states at certain 
rational fractional values of Landau level filling factor 
v = nh/eB. Experimentally, the "magnetic field induced 
electron solid" (MIES) 3 forms at sufficiently small u, fol- 
lowing the termination of FQH states at low T. It is 
deformed and pinned by disorder in the system, thus an 
insulator (in DC transport), and has a characteristic res- 
onance (see, for example, Ref. 14) in its frequency (/) 
dependent, real diagonal conductivity (Re[<7 xx (/)]) due 
to the "pinning mode" of domains of the elastic solid 
oscillating collectively around the disorde»i^iiS*ii. 

The melting of such a MIES has been studied in 
various experiments 18 - 19 ' 20 ' 21,22 - 23 and it was commonly 
presume d 8 i 18 ' 20 i 22 that at a fixed v the melting should be 
similar to that of a classical electron solid (the expected 
exact ground state for a 2DES at infinite B) and thus 
T m would be determined by n, as T cm is. In this letter, 
we show unambiguously that this is not the case, that 
is, the melting of a MIES does not behave as a classical 
electron solid, and in fact in any given sample (whose 
intrinsic disorder is fixed), T m is determined by v, not n, 
and is unrelated to T cm . 

In our experiments we have studied the In- 
dependence of the pinning mode resonance of the MIES 
in its R.e{o~ xx (f)] spectrum measured by microwave 
spectroscopy 2 ^. No resonance is observed when T is 
raised above some characteristic T m , taken as the melt- 
ing T of the electron solid. By systematically measuring 
T m while varying both n and B, we found that within 
the experimental resolution, T m in a given sample is only 
a function of v (i.e., T m (n,B) = T m (n/B)) down to v 
as small as ~0.05 attained in our experiments. Although 
T m is generally sensitive to n at fixed B, T m is insensitive 
to n at fixed v. Since v — nh/eB — 2(Zs/r) 2 , where the 
magnetic length Ib = y/H/eB is a measure of the size of 
the single electron wavefunction and r = 1/^/irn is the 
mean separation between the electrons, our findings re- 
flect the quantum nature of the 2D electron solid formed 
at finite B and demonstrate that its melting (T m ) is ef- 
fectively controlled by the inter-electron quantum corre- 
lation, which depends on Ib/t. 

We have studied two different 2DES samples. Sam- 
ple 1 is a GaAs/AlGaAs heterojunction. Sample 2 is a 
15nm-wide AlGaAs/GaAs/AlGaAs quantum well (QW). 



By backgating and/or different cooldowns, the electron 
densities (n) in both samples can be tuned to various ex- 
tents (to be specified below). Typically at their respec- 
tive as-cooled n, sample 1 has mobility /i~6x f 6 cm 2 /Vs 
, sample 2 has ^lxlO 6 cm 2 /Vs . 

The inset of Fig. IHb) shows a schematic local cross 
section of a typical sample (not to scale). A metal film 
coplanar waveguide (CPW), lithographically deposited 
on the surface, enables the measurement of Re[cr ra (/)] 
of the 2DES. A network analyzer generates a microwave 
signal propagating along the CPW and coupling capac- 
itively to the 2DES (located some 0.1-0.5 fim below the 
surface) . The relative power absorption (P) by the 2DES 
is measured and can be related 2 - 4 to Refer^] of the 2DES 
as P = ex.p((2l Z Q / w)Re[a xx \) , where I and w are the to- 
tal length and slot width of the CPW respectively and 
Zq=50£1 is the CPW characteristic impedance. CPW of 
meander shapes^ are commonly used to obtain larger 
geometric ratios (21 /w), and therefore to increase the 
strength of the signal (P). The sample is mounted on 
a metal block that is kept in good thermal equilibrium 
with the mixing chamber of a dilution refrigerator during 
the measurements. The microwave input is kept in the 
low power limit, by reducing the power till the signal (P) 
no longer changes. A negative voltage between a back- 
gate (located ^200 /am below the surface) and the 2DES 
enables in-situ reduction of n from the as-cooled values. 

Fig.^shows T-dependence of the microwave resonance 
of the electron solid in sample 1 and the determination 
of T m at two different values of n with u fixed at a rep- 
resentative value of 0.128. In Fig.[TJi, n=5.6xl0 10 cm" 2 
(5=18 T) and the Re[a xx (f)] spectrum displays a clear 
resonance with peak frequency (/ p k) near 600 MHz at low 
T(~50 mK). As T is increased, the resonance weakens. 
Its / p k at elevated T also decreases slightly from the base 
T (50 mK) value, but by no more than 20%, indicating 
that there is no significant loss of pinning of the elec- 
tron solid due to the increase of the temperature. The 
resonance disappears into noise background at ~250 mK, 
which we take as the melting T (T m ) of the electron solid. 
The inset of Fig. QJi shows that the resonance amplitude 
(obtained from a Lorentzian fit) extrapolates to zero at 
a similar T m . In Fig.^ 3 , n has been reduced to 2.1 xlO 10 
cm" 2 (using a backgate voltage of —300 V) while B is 
also reduced to 6.86 T to keep v the same value as in 
Fig. ^,. The lower n reduces the electron-electron inter- 
action relative to the electron-disorder interaction, and 
the pinning resonance of the electron solid now occurs (at 
low T) near 1.6 GHz, which is significantly higher than 
that in Fig.^i. However, the resonance disappears above 
a T m ~ 270 mK, which is similar to the T m in Fig. QJi, 
within the experimental uncertainty in determining T m 
(typically -10%). 

We have measured T m in 4 different cooldowns of sam- 
ple 1 for many combinations of n and B, at which a reso- 
nance from the MIES can be detected. We plot all these 
T m data as a function of v in Fig. [21 Sample 1 has a par- 
ticularly large range of tunable n, from ~1.2— 8.1xl0 10 
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FIG. 1: Temperature (T) dependence of the microwave spec- 
tra of the 2D electron solid measured from sample 1 at two 
different densities (n) with the same Landau filling factor 
(v=nh/eB). a, T-dependence of the spectra at n=5.6xl0 10 
cm -2 with */=0.128 (magnetic field B=18 Tesla). Spectra at 
a series of representative T's are shown and offset for clar- 
ity. The pinning resonance of the electron solid observed at 
T=50 mK is seen to weaken with increasing T and disappear 
at ~250 mK, taken as the melting temperature (T m ) of the 
electron solid. Inset shows the amplitude of the resonance ex- 
trapolates to zero at the similar T m . b, T-dependence of the 
spectra at n=2.1xl0 10 cm" 2 with i/=0.128 (B=6.86 Tesla). 
The low T resonance disappears at a T m similar to that shown 
in a (n=5.6x 10 10 cm" 2 ), despite here that n has been reduced 
by more than a factor of 2 (using a negative voltage between 
a backgate and the 2DES). Inset shows a schematic of the 
sample. The dark regions on the top surface represent the 
coplanar wave guide (CPW). The backgate allows the in-situ 
change of n. 



cm" 2 , covering a v range from ~0.21 down to ~0.03. 
Generally, we found T m to be sensitive to n or B sepa- 
rately (when fixing one and changing the other). On the 
other hand, near similar v (changing both n and B while 
fixing n/B), we have always found T m to be insensitive 
to n (or B), within the experimental errors in T m . Dif- 
ferent cooldowns can vary T m by up to ~15% (at similar 
v) but this does not affect our conclusion. We thus find 
T m to be mainly determined by v, so that T m vs v as 
plotted in Fig. [21 defines a melting curve for the electron 
solid in sample 1. A linear fit of T m vs v gives a guide 
to the eye shown as the dashed line, which lies within 
20% from all the (T m , v) data points and within 10% 
from a majority (70%) of them. In the inset of Fig. [21 
we plot the "reduced"^*^ t m = T m /T cm from two sim- 
ilar cooldowns versus v, where T cm = e 2 ^/n/(A r KtkB^) 
is the melting T of a classical 2D electron solid defined 
earlier (with the value 78 T=127). In contrast to T m , 
t m can vary significantly (sometimes by a factor of 3) 
at similar v. Thus t m versus v does not give a well 
defined melting curve for the electron solid, confirming 
that T m is not determined by n or T cm . We have also 
checked that neither T m nor t m plotted against either n 
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FIG. 2: (T m ,v) phase diagram for the electron solid in sam- 
ple 1. The T m 's are measured in a total of 4 cooldowns 
(shown as diamonds, circles, squares and triangles), over a 
wide range of densities (n—1.2— 8.1xl0 10 cm" 2 ) and mag- 
netic fields. Within the experimental uncertainty, T m versus v 
gives rise to a well-defined melting curve of the electron solid. 
The dashed line is a guide to the eye, obtained by a linear fit 
through all the data. Typical error bars in T m are less than 
10%. The inset shows the "reduced" t m versus v from two 
cooldowns. tm is defined as T m normalized by the classical 
2D electron solid melting temperature T cm — e 2 y/n/(AitekBT) 
(where we take F=127). t m versus v does not result in a well- 
defined melting curve, indicating that the meting temperature 
T m is not determined by n or T cm . 



or B gives any well defined melting curve. We further 
notice that our measured T m is also unrelated to the 
hypothetical "thermal depinning" temperature (Tdcpin) 
of domainsi5ii£ii£ of the electron crystal. In fact, one 
can estimate Tdepin^(l/fcs)m-eWo^ 2 (^/ a ) 2(xn3 ^ 2 (where 



=2nf pk eB/j 



3 / 2 characterizes the "effective pin- 



ning potential"— (we note that wo is also much larger 
than ksTm/H), £ is the disorder correlation length, Loan 
is the Larkin domain sizei 5 ^i§*i 7 » and a is the electron 
crystal lattice constant), in contrast to T m , which is de- 
termined by v, and as seen in Fig. [21 typically decreases 
with increasing n if B is fixed. 

Fig. |31 shows the (T m ,v) melting curve measured on 
sample 2 (15nm-wide QW). Sample 2 has a tunable 
rt=2.7-4xl0 10 cm" 2 . Likely due to the relatively narrow 
confinement of the 2DES in the QW, sample 2 enters the 
solid phase for ^<0.3 and the typical / p k observed is ^6-8 
GHz. In this sample we found T m again to be controlled 
by v, though the value of T m is higher than in sample 1 
at similar v. 

In each of the two samples, at a fixed high-S, T m typi- 
cally increases with decreasing n. This is opposite to the 
classical behavior, where T cm decreases with decreasing 
n. At a fixed v, our observed T m is insensitive to n, and 
T m in each sample is a well defined function of v. Re- 
cent theories 2 ^* 2 ^ have suggested that many-body quan- 
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FIG. 3: (T m ,t/) melting curve of the electron solid in sample 
2, a narrow QW of width 15 nm. Sample 2 has a tunable 



n=2.7-4xl0 1 
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cm and enters an electron solid phase at u< 



turn correlation between electrons can still be important 
in the solid phase terminating the FQH states at high B. 
The parameter that captures such quantum correlation 
is the Laudau filling factor z/=2(^/r) 2 , which also de- 
termines the (single-particle) wavefunction overlap^ 7 , be- 
tween neighboring electrons. This overlap (I s ), given 27 
by e -( 2 /v3)7r/^ i s j n f ac t quite small in the high-i?, low- 
v regime of the MIES. Our data show that T m is quite 
strongly dependent on v (for example, T m changes by 
almost a factor of 2 when v is changed by a factor of 
2). Therefore the dependence of T m on v is unlikely 
to result from a T m that is determined by the Coulomb 
Hartree energy Eh (which is showr*2i to have very weak 
^-dependence in this low-i/ regime, where the wavefunc- 
tion overlap I s has negligible contribution to Eh)- Our 
findings indicate that the melting of MIES is effectively 
controlled by the inter-electron quantum correlation^!, 
through v. The well-defined (T m ,v) melting curve we 
obtained constitutes the phase boundary between a quan- 
tum solid and a correlated quantum liquid 15 in each sam- 
ple. 

The T m we measured in both samples are of 
similar order of magnitude to those in other sam- 
ples measured previously with various experimental 
techniquea 18 i 19 ' 20 i 21 i 22 i 23 . We have noticed that at simi- 
lar v, sample 2 (narrow QW) has higher T m than sample 
1 (heterojunction). It has been suggested^ that the rel- 
evant disorder in semiconductor samples that pins the 
MIES comes mostly from the interfaces vertically confin- 
ing the 2DES (this is consistent with the relatively high 
/ p k observed in the narrow QW sample). Our findings 
thus suggest that a 2D electron solid subject to stronger 
pinning disorder melts at higher T m , an effect that has 
been predicted earlier 28 . We notice, interestingly, that 
studies on the effects of pinning disorder or geometric 



confinement on T m in other solids usually find the op- 
posite behavior. For example, both the vortex solid in 
a high-temperature superconductor subjected to artifi- 
cial pinning centers^ and a helium solid in a porous 
glass^i show a depression of T m with the added disor- 
der. Disorder is generally unavoidable in semiconductor 
samples, and as seen here can strongly influence the val- 
ues of T m measured, which may not be those of an ideal, 
clean Wigner crystal. On the other hand, the observa- 
tion that in any given sample (thus with fixed intrinsic 
disorder), T m is not sensitive to n (nor to B) if v=nh/eB 
is fixed (but is sensitive to n if B is fixed) is difficult to 
explain in a classical picture of, for example, interplay 
of disorder and n (screening). The fact that v being the 
controlling variable for T m is found in samples of quite 
different disorder (including sample 1, in which one can 
be infer from the observed pinning resonance that the 
Wigner crystal there can possess substantial correlation 
lcngt h 14 i 15 i 16 i 17 ) suggests that such a behavior (T m be- 
ing controlled by v) arises from some mechanism in which 
many-body quantum correlation between electrons may 
be important. The exact role of disorder (perhaps involv- 
ing its interplay with quantum correlation) in determin- 
ing the value and behavior of T m remains to be better 
understood. 

It is also interesting to compare the quantum nature 
of our 2D electron solid to the quantum solids of he- 
lium. In a helium solid, the size of atoms is fixed 
by nature and the quantum parameter is the De Boer 
parameteni A~h/(a\/Mv) (where M is the atomic mass, 
a the inter-atomic distance and v the inter-atomic po- 



tential strength), which is fixed at fixed n. The T m of a 
helium solid is only determined by n. In the case of a 2D 
electron solid formed in high B, the size of single electron 
wavefunction (Ib) is readily tunable by B, independently 
of n. The quantum parameter here is v=nh/eB and we 
have found that T m of such an electron solid in a given 
sample (with fixed disorder) is controlled by v rather than 
n. 

Acknowledgement Financial support of this work 
was provided by DOE grant no. DE-FG02- 
05ER46212, the NHMFL in-house research pro- 
gram, and AFOSR. The spectroscopy measure- 
ments were performed at the National High Mag- 
netic Field Laboratory, which is supported by NSF 
Cooperative Agreement No. DMR-0084173 and by 
the State of Florida. We thank Glover Jones, Tim 
Murphy and Eric Palm at NHMFL for experimen- 
tal assistance. We also thank S. T. Chui, R. L. 
Willett, Kun Yang and C. C. Yu for inspirational 
discussions, and thank H. A. Fertig and P. B. Lit- 
tlewood in particular for discussions regarding the 
thermal depinning temperature of the Wigner crys- 
tal. 

Competing Interests The authors declare that they 
have no competing financial interests. 

Correspondence Correspondence should be addressed 
to Y.P.C. (email: yongchen@rice.edu). 



* Current address: Richard Smalley Institute for Nanoscale 
Science and Technology and Dept. of Physics, Rice Uni- 
versity, Houston TX 77005, USA 

' Current address: Dept. of Physics, University of Maryland, 
College Park MD 20742, USA 

* Current address: School of Electrical and Computer En- 
gineering, Purdue University, West Lafayette, IN 47907, 
USA 

1 Wilks, J. The Properties of Liquid and Solid Helium (Ox- 
ford, 1967). 

2 Kim, E. & Chan, M. H. W. Probable observation of a 
supersolid helium phase. Nature (London) 427, 225-227 
(2004). 

3 Shayegan, M. Case for the magnetic-field induced two- 
dimensional Wigner crystal. In Das Sarma, S. & Pinczuk, 
A. (eds.) Perspectives in Quantum Hall Effects, chap. 9 
(Wiley and Sons, New York, 1997). 

4 Grimes, C. C. & Adams, G. Evidence for a liquid-solid 
phase transition in a classical, two-dimensional sheet of 
electrons. Phys. Rev. Lett. 42, 795-798 (1979). 

5 Wigner, E. P. On the interaction of electrons in metals. 
Phys. Rev. 46, 1002-1011 (1934). 

6 Thouless, D. Melting of the two-dimensional Wigner lat- 
tice. J. Phys. C 11, L189-L190 (1978). 

7 Morf, R. H. Temperature dependence of the shear mod- 
ulus and melting of the two-dimensional electron solid. 



Phys. Rev. Lett. 43, 931-935 (1979). 

8 Chui, S. T. & Esfarjani, K. Finite-temperature two- 
dimensional Wigner transition. Phys. Rev. B. 44, 11498- 
11501 (1991). 

9 Lozovik, Y. E. & Yudson, V. I. Crystallization of a two- 
dimensional electron gas in a magnetic field. JETP Lett. 
22, 11-12 (1975). 

10 Fukuyama, H., Platzman, P. M. & Anderson, P. W. 
Two dimensional electron gas in a strong magnetic field. 
Phys. Rev. B. 19, 5211-5217 (1979). 

Laughlin, R. B. Anomalous quantum Hall effect: An in- 
compressible quantum fluid with fractionally charged ex- 
citations. Phys. Rev. Lett. 50, 1395-1398 (1983). 

12 Kivelson, S., Kallin, C, Arovas, D. P. & Schrieffer, J. R. 
Cooperative ring exchange theory and the fractional quan- 
tum Hall effect. Phys. Rev. B. 36, 1620-1646 (1987). 

13 Tsui, D. C, Stormer, H. L. & Gossard, A. C. Two- 
dimensional magnetotransport in the extreme quantum 
limit. Phys. Rev. Lett. 48, 1559-1562 (1982). 

14 Ye, P. D. et al. Correlation lengths of the Wigner-crystal 
order in a two-dimensional electron system at high mag- 
netic fields. Phys. Rev. Lett. 89, 176802 (2002). 

15 Chitra, R., Giamarchi, T. & Le Doussal, P., Dynamical 
properties of the pinned Wigner crystal, Phys. Rev. Lett. 
80, 003827 (1998); Pinned Wigner crystals, Phys. Rev. B 
65, 035312 (2002). 



16 Fertig, H. A. Electromagnetic response of a pinned Wigner 
crystal. Phys. Rev. B. 59, 2120-2141 (1999). 

17 Fogler, M. M. & Huse, D. A. Dynamical response of a 
pinned two-dimensional Wigner crystal. Phys. Rev. B. 62, 
7553-7570 (2000). 

18 Glattli, D. C. et al. Experiments on ordering in classical 
and quantum 2d electron systems. Surf. Set. 229, 344-351 
(1990). 

19 Goldman, V. J., Santos, M., Shayegan, M. & Cunning- 
ham, J. E. Evidence for two-dimensional quantum Wigner 
crystal. Phys. Rev. Lett. 65, 2189-2192 (1990). 

20 Williams, F. I. B. et al. Conduction threshold and 
pinning frequency of magnetically induced Wigner solid. 
Phys. Rev. Lett. 66, 3285-3288 (1991). 

21 Paalanen, M. A. et al. Electrical conductivity and Wigner 
crystallization. Phys. Rev. B. 45, 13784-13787 (1992). 

22 Goldys, E. M. et al. Magneto-optical probe of two- 
dimensional electron liquid and solid phases. Phys. Rev. B. 
46, 7957-7960 (1992). 

23 Kukushkin, I. V. et al. Wigner solid vs. incompressible 
Laughlin liquid: phase diagram derived from time-resolved 
photoluminescence. Europhys. Lett. 23, 211-216 (1993); 
Kukushkin, I. V. et al. Evidence of the triangular lattice 
of crystallized electrons from time resolved luminescence, 
Phys. Rev. Lett, 72, 3594 (1994). 

24 Engel, L. W., Shahar, D., g. Kurdak & Tsui, D. C. Mi- 



crowave frequency dependence of integer quantum Hall ef- 
fect: Evidence for finite-frequency scaling. Phys. Rev. Lett. 
71, 2638-2641 (1993). 

25 Yi, H. & Fertig, H. A. Laughlin- Jastrow-correlated Wigner 
crystal in a strong magnetic field. Phys. Rev. B. 58, 4019- 
4027(1998). 

26 Chang, C. C, Jeon, G. S. & Jain, J. K. Microscopic ver- 
ification of topological electron-vortex binding in the low- 
est landau-level crystal state. Phys. Rev. Lett. 94, 016809 
(2005). 

27 Maki, K. & Zotos, X. Static and dynamic properties of a 
two-dimensional Wigner crystal in a strong magnetic field. 
Phys. Rev. B. 28, 4349-4356 (1983). 

28 Tsukada, M. Two-dimensional crystallization of the elec- 
trons in MOS structures induced by strong magnetic field. 
J. Phys. Soc. Jpn. 42, 391-398 (1977). 

29 Christensen, H. O. Confinement effects on freezing and 
melting. J. Phys. C 13, R95-R133 (2001). 

30 Paulius, L. M. et al. Evolution of the vortex phase 
diagram in YBa2Cu307-<s with random point disorder. 
Phys. Rev. B. 61, R11910-R11913 (2000). 

31 Beamish, J. R., Hikata, A., Tell, L. & Elbaum, C. So- 
lidification and superfluidity of 4 He in porous vycor glass. 
Phys. Rev. Lett. 50, 425-428 (1983). 



